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Abstract 

We study the pion production in heavy-ion collisions at SIS energies within 
the coupled resonance BUU model. The standard calculation overpredicts the 
pion multiplicity in the central Au+Au collision at 1 A GeV by about a factor 
of 2, which is interpreted as a signature for the quenching of the resonance 
production/absorption in nucleon-nucleon collisions at high nuclear density. 
Results of the calculations taking into account a specific form of quenching 
are presented for pion spectra from light and heavy systems as well as for the 
pion in-plane and out-of-plane flow. 

I. INTRODUCTION 

Modifications of hadrons in dense nuclear matter are currently in the center of theoretical 
studies 0,0] ■ 111 particular, the multiple pion production in heavy-ion collisions offers an 
interesting possibility to look at the properties of the baryon resonances in nuclear matter 
PHlT|. It is well established (c.f. ||[| ) that at beam energies of 1 4- 2 A GeV pions are 
produced mostly in the two-step process NN — > NR, R — > Ntt, where R stands for the 
A(1232) or higher baryon resonances. However, there is an overall tendency of transport 
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models explicitly propagating resonances to overpredict the pion multiplicity in heavy-ion 
collisions at the beam energy ~ 1 A GeV (c.f. ||12|| ). Furthermore, taking into account 
in-medium pion dispersion relations even increases the pion yield |13|-|i"5 . 



In the present work we study the effect of possible in-medium modifications of the reac- 
tion rates 

NN <-> NR (1) 

on pion production. A description of the experimental data on the pion multiplicities re- 
quires an in-medium reduction of the rates (|]). We propose, therefore, a density-dependent 
modification of the transition rates ([!]), which gives a reasonable description of the pion 
multiplicities in the systems C+C, Ni+Ni and Au+Au at the beam energy of l-i-2 A GeV. 
We will show, moreover, that this modification is consistent with other pion observables, 
i.e. transverse mass and kinetic energy spectra, angular distributions and in-plane collective 
flow. 

The structure of the paper is as follows: In Sect. II we discuss the effect of the resonance 
production/absorption quenching on the pion production. Sect. Ill contains the results on 
the collective pion flows, while in Sect. IV the summary and conclusions are given. 

II. DENSITY DEPENDENT TRANSITION RATES NN <-» NR 

In our numerical calculations we have used the BUU model in the version described in 
Ref. [[UJ. The nucleons, A(1232) and higher (up to D 35 (2350)) baryon resonances, A and E 



hyperons, mesons 7r, t], u, K, p, a are explicitly propagated in this BUU implementation. 
The model has been successfully applied to the calculation of nucleon collective flow observ- 



ables in heavy-ion collisions at energies from ~ 100 AMeV up to ~ 2 A GeV | I7[| . A good 
description of " more elementary" 7 and e~ induced reactions has been also reached within 
the model (see Refs. |B|,@). 

In this work we concentrate on pion production in nucleus-nucleus collisions continuing 
the studies initiated in Refs. |6,[lTl. All calculations in the present paper are performed using 



the soft momentum-dependent mean field (SM) with incompressibility K = 220 MeV, which 



gives the best overall description of the nucleon collective flow variables [[[7]]. The BUU set 
of parameters, which includes the SM mean field as well as vacuum resonance parameters 
and vacuum elastic and inelastic cross sections, will be denoted as "standard" below. The 
modifications of the standard set will be discussed in the text explicitly. 

First, we compare the calculated pion multiplicities in central Au+Au collisions at 1 
A GeV with the experimental data. In Fig. 1 the dashed line shows the pion number vs 
time in the case of the standard calculation. The pion multiplicity asymptotically reaches a 
value of ~ 60 which is almost a factor of 2 larger than the experimental value of ~ 34 (c.f. 
Ref. fL2|). Pions are produced dominantly by the baryon resonance decays. Therefore, the 
pion multiplicity is very sensitive to the number of resonances present in the system. The 
differential production cross section of a resonance R with mass M in a nucleon-nucleon 
collision, averaged over the spins of initial particles and summed over the spins of final 
particles, is given as 

d -^^ = ^^AM*umj^), (2) 



where \A4\ 2 is a spin-averaged matrix element squared, P12 and P34 are cm. momenta of 
incoming and outgoing particles, e is the total cm. energy, A(M 2 ) is the spectral function 
of the resonance, and Jr is the spin of the resonance. The inverse reaction cross section is 

da N3Ri ^ NlN2 _ 1 1-7712 Pjj v 4 

where, due to detailed balance, the same matrix element appears as in Eq. (fj), C\i = 2 if 
the nucleons N\ and N2 are identical and G\i = 1 otherwise. 

The matrix element which appears in Eqs. (@),© should be evaluated in nuclear matter. 
Instead of using any model dependent parametrizations for the in- medium matrix element, 
we adopt a simple relation 



\M\ 2 = K(p)\M vac \ 2 , k(0) = 1, (4) 
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where |-M mc | 2 is the vacuum matrix element squared and n(p) is some density-dependent 
function to be determined by comparison with the experimental data. The in-medium 
matrix elements for the non-resonant (background) one-pion production/absorption N\N 2 <-> 
double A(1232) production/absorption N\N 2 A 3 A 4 and resonance scattering on 
nucleons N\R 2 — > N 3 R4 are also related to their vacuum values by Eq. (Q) with the same 
universal function n(p). 

We address now the pion multiplicity in central Au+Au collisions at 1.06 A GeV using 
various choices of the function n(p). As a first trial, we simulate an enhancement of resonance 
production in NN collisions [0] assuming a linear dependence on density as in Refs. |TT| , |20 



«(» = 1 + /?-, (3>0, (5) 
Po 

where po — 0.16 fm -3 is the equilibrium nuclear matter density. We see from Fig. 1 
(dotted line), that the amplified matrix element with /t(p) = 1 + 3p/po leads to an about 
20% reduction in the final pion multiplicity which, however, is still above the experimental 
abundance. This decrease originates from the amplified - by detailed balance - absorption 
of resonances in the expansion stage of the system. A further increase of the parameter f3 
practically does not reduce the pion number any further (see Fig. 3 for (3 > 0). 

An essential reduction of the pion multiplicity can only be reached by a quenching of the 
resonance production/absorption processes in nuclear matter at high density. To this aim 
we apply a piece-wise linear function of the form 

k( P ) = min(l,max(0,l + /3(p/p -l))) , P<0, (6) 

which is shown by the solid line in Fig. 2 for the particular case (3 = —2. Within the 
parametrization (||), the in-medium effects vanish at subnuclear densities (k = 1), but 
above po the nuclear medium creates the shadowing of binary collisions involving resonances 
leading, effectively, to a quenching of these collisions for p > (1 — l/j3)p . In Fig. 3 (see the 
part of the BUU curve for (3 < 0) we report the number of produced negative pions vs (3 in the 
case of the parametrization (§). The proper pion multiplicity is reproduced for (3 = —2; we 



will denote this calculation as "quenched" in the following. For this calculation, resonances 
do not experience any elastic or inelastic scatterings in high density nuclear matter. They 
can, however, decay or be produced in the processes R <->• Ntt. In the quenching scenario, the 
resonance production in NN collisions happens only at relatively low density p < 1.5/?o> i- e - 
dominantly in the overlapping tails of the density profiles of colliding nuclei. This reduces 
the maximum number of A-resonances to about 35% as compared to a standard calculation 
(c.f. solid and dashed lines in the lower panel of Fig. 1). 

Before discussing other pion-related observables we check, how the resonance quenching 
influences the underlying nucleon dynamics. Fig. 4 shows the proton rapidity distribution 
and the transverse momentum of protons in the reaction plane as a function of rapidity. 
Due to a lower nucleon inelasticity within the quenching scenario, the nuclei show less 
stopping, i.e. the rapidity distribution becomes somewhat wider and, since also the pressure 
is decreased, a slightly smaller transverse in-plane flow is generated making the agreement 
of the standard SM calculation JET] with the proton flow data worse (see also Fig. 11). This 
problem, however, can be resolved by either increasing the incompressibility or the stiffness 
of the momentum-dependent interaction. The latter modifications influence only weakly the 
pion multiplicity and, therefore, are not studied in this work. 

In Fig. 5 we present the pion multiplicity as a function of the participant number A part for 
Au+Au collisions at 1.06 A GeV in comparison with the FOPI and TAPS data collected in 



Ref . [0 . We have determined the participant number for a given impact parameter using the 
participant-spectator model with sharp nuclear density profiles. The quenched calculation is 
in a good agreement with the data for central collisions; in peripheral collisions the medium 
effects are smaller, the difference between the standard and the quenched calculation vanishes 
gradually and, thus, the multiplicity is still somewhat overpredicted. This discrepancy in 
peripheral collisions could arise from different definitions of the participant number A part in 
our calculations and in Ref. ||12| , since the impact parameter can not be measured directly 
in the experiment. 

Fig. 6 shows the beam energy dependence of the charged pion multiplicity for central 



Ni+Ni collisions compared with the FOPI data from Ref. |2T|| . The quenched calculation 
describes the data rather well. With increasing beam energy the difference between stan- 
dard and quenched calculations reduces, since the nonresonant string mechanism for pion 
production becomes dominant, which is the same in both calculations. We recall that in 
our transport approach strings are excited in NN collisions for invariant collision energies 
y/s > 2.6 GeV corresponding to the beam energy of 1.73 A GeV if the Fermi motion is 
neglected. 

Fig. 7 shows the cm. kinetic energy spectra of charged pions at c . m . = 130° for Au+Au 
at 1.06 A GeV. The quenching improves the agreement with the data somewhat producing 
a more concave shape of the spectrum for E%-£' > 0.1 GeV. However, the calculated cross 
sections in this region are still too large, while at low kinetic energy the experimental data 
are underpredicted by the quenched calculation. 

A similar tendency is visible in the case of pt-spectra at midrapidity shown in Fig. 8: 
the improvement by the quenching mechanism at high p t is evident, but at low p t too few 
pions are produced. This leads also to the underprediction of the midrapidity pion yield by 
the quenched calculation. Therefore, since the total pion multiplicity agrees with the data 
(or slightly overpredicts for peripheral collisions, see Fig. 5), the quenched calculation will 
produce some excess of pions at projectile/ target rapidities. We will, however, discard here 
the comparison of our results with the FOPI data on the pion rapidity distribution from 



Ref. | 12|| . This comparison would require the treatment of our BUU-events in exactly the 
same way as the data, which is out of scope of this work. 

The reason for the underprediction of the low-p t pion yields by the quenched calculation 
is, probably, an overdamping of the low mass resonance production in NN collisions at high 
density by the quenching factor (|6]). Indeed, a smaller medium effect is expected for higher 
kinetic energies of the final particles in analogy to Pauli blocking p3| . In this exploratory 
work we are not aiming at a detailed description of the pion energy spectra, but, rather, 
want to demonstrate the gross effect of the density modified matrix element on various 
pion observables. It is straightforward, however, to introduce a \fs and resonance-mass- 
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dependent modification factor k in Eq.(^). Another physical effect is the off-shellness of the 
pion in the nuclear medium due to large nN — > A cross section, which leads to a broad pion 
spectral function. Thus a dynamical evolution of the hadron spectral functions in BUU will 
modify the low momentum spectra sensitively ||24j| . 



In Fig. 9 we present the cm. polar angle distributions of charged pions. The shape of 
these spectra is similar to the one for free-space pion production in NN collisions. There is 
a good overall agreement of the quenched calculation with the data. 

For the collisions C+C at 0.8-=-2 A GeV (Fig. 10) the calculated mt-spectra of 7r°'s are 
practically identical in the standard and quenched calculations. This is expected, since the 
in-medium modifications are smaller in the lighter system. 

III. COLLECTIVE FLOW 

In earlier studies (c.f. [f4|,p6f ) , the pion antiflow and squeeze-out have been explained 
by the rescattering of pions on the spectator nucleons, which leads to the shadowing phe- 
nomenon. 

Fig. 11 shows the proton, 7r + and ir~ in- plane flow (:= d < p x > /dY<® at Y<® = 0) 
for the Au+Au system at 1.15 A GeV as a function of impact parameter in comparison 



to the data p7|| . The best overall agreement with the data is reached within the standard 
SM calculation. For ir ± , the quenched and standard results are almost indistinguishable 
within statistics. However, the proton flow is underpredicted by the quenched calculation 
(see discussion in the previous section). The difference in flows for 7r + and tt~ is due to the 
Coulomb interaction with protons which increases (decreases) the antiflow for 7r + (7r~) for 
large impact parameter. 

In Fig. 12 we show the squeeze-out ratio R N := {N w+ {90°) + A^ + (-90°))/(A^+(0°) + 
^+(180°)) as a function of the transverse momentum for positive pions at midrapidity in 



comparison to the data [28]. The squeeze-out ratio grows with p t since faster moving pions 



in transverse direction reach the spectators earlier and rescatter in analogy to the case of 
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the nucleon squeeze-out fl7fl . A peculiar feature of the -Rjv(pt)-dependence for positive pions 
is its saturation at relatively low p t ~ 0.3 GeV/c, while the nucleon squeeze-out ratios 



increase monotonically with p t (c.f. [EMROf). This behaviour of the pion squeeze-out ratio is 



related to the resonance nature of ttN scattering ||26|| . Our calculations, both standard and 
quenched, underpredict the absolute value of R^, but the shape of the function R N (p t ) is in 
a qualitative agreement with the data. The reason for this discrepancy lies, most probably, 
in the absence of any momentum dependence in the pion potential. 

IV. SUMMARY AND CONCLUSIONS 



Within the transport BUU model |16j we have studied the influence of the in-medium 
modifications of the reaction rates (|]) on pion production in heavy-ion collisions around 
1 A GeV. The pion multiplicity in the system Au+Au at 1 A GeV can be consistently 
described assuming a quenching of the resonance production/absorption in NN collisions at 
high nuclear densities. The quenching generally improves the description of standard BUU 
calculations with experiment on other pion observables: pt-, m t - and kinetic energy spectra, 
in-plane and out-of-plane collective flows. 

The physical reason for the in-medium reduction of the resonance production rates, that 
seems to be needed to describe the pion data, remains open. This is in particular so because 
there are theoretical predictions that the A excitation cross section should increase with 
density ||19||. A collisional broadening of the A-resonance, that was not taken into account 



in Ref. |19|], could, in principle, provide the obviously needed reduction, but the calculations 



of Ref. |3T| , that were based on reasonable collision rates, have shown that this broadening 
is too small. All these models have in common that they work with free coupling constants 
and thus neglect any in-medium vertex corrections. Good calculations of inelastic nucleon- 
nucleon scattering in nuclear matter would be needed to clarify this problem. 
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FIGURE CAPTIONS 



Fig. 1 Pion number (upper panel) and A-resonance number (lower panel) vs time for a 
central collision Au+Au at 1.06 A GeV. Standard calculations are represented by 
dashed lines. The dash-dotted and dotted lines show the cases of in-medium amplified 
reaction rates (JTp for (3 — 1 and (3 — 3, respectively, while the solid lines correspond to 
the quenched calculation ((3 = —2). The experimental data from Ref. \T2\ are within 
the dashed band in the upper panel. 

Fig. 2 Density dependence of the modification factor n(p) of the matrix element for reso- 
nance production/absorption (see Eq.(fD). The notation of lines is the same as in Fig. 
1. 

Fig. 3 Number of negative pions for the central collision Au+Au at 1.06 A GeV as a 
function of the coefficient (3 in the density dependence of the factor n (see Eqs.(|]),(|D). 
The part at /3 < corresponds to the parametrization (|6|), while n(p) = 1 + Pp/po for 
(3 > 0. The dashed band represents the experimental data from Ref. [12| . 



Fig. 4 Proton rapidity distribution (upper panel) and rapidity dependence of the proton 
transverse momentum in the reaction plane (lower panel) for Au+Au collision at 1.06 
A GeV and impact parameter b=6 fm. The rapidity is normalized to the projectile 
rapidity in the center-of-mass system: = (y/y P roj)c.m.- 

Fig. 5 7r~, 7T° and it + multiplicities (left, central and right panels, respectively) vs partici- 
pant number A part for the system Au+Au at 1.06 A GeV. The experimental data are 



from Ref. 12 



Fig. 6 7i and 7r + multiplicities (left and right panels, respectively) in central collisions of 
Ni+Ni as a function of the beam energy. The BUU calculations are performed for b=l 



fm. The experimental data are from Ref. [21 . 
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Fig. 7 Cm. kinetic energy inclusive spectra of negative and positive pions (upper and lower 
panels, respectively) at c . m . = 130° for Au+Au at 1.06 A GeV. The experimental data 



are from Ref. 12 



Fig. 8 7r~, 7T° and 7r + transverse momentum inclusive spectra at midrapidity for Au+Au 
at 1.06 A GeV. The spectra are calculated in the rapidity intervals = —0.2 + 0.2 
for 7r ± and = —0.25 + 0.21 for 7r°, but normalized to a rapidity interval dy = 1 



following Refs. [12,22]. The experimental data are from Refs. [p 



Fig. 9 Cm. polar angle distributions of charged pions in Au+Au collisions at 1.06 A GeV. 
Left and right panels show, respectively, the ir~ and tt + distributions for the cm. 
kinetic energy Ef^' > 40 MeV, while the central panel presents the ir~ distribution 
for > 365 MeV. The experimental data are from Ref. [0 . 

Fig. 10 Transverse mass spectra of 7r°'s in selected rapidity intervals near midrapidity 
(following Ref. |2§) for C+C collisions at 0.8, 1 and 2 A GeV: y lab = 0.42 -=-0.74, 0.42 + 



0.74 and 0.8 + 1.08, respectively. The thick solid lines show exponential fits to the 
experimental data from Ref. P 



Fig. 11 In-plane flow as a function of impact parameter for protons (upper panel), 7r +, s 
(middle panel) and 7r~'s (lower panel) in comparison with the experimental data from 
Ref. |3 for Au+Au at 1.15 A GeV. 



Fig. 12 Transverse momentum dependence of the squeeze-out ratio Rn for positive pions 
for semicentral Au+Au collisions at 1 A GeV. The calculated results are impact- 
parameter-weighted in the range 6 = 5 + 7 fm. The experimental data are from Ref. 
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